The Astrophysical Journal, 781:107 (llpp), 2014 February 1 

© 2014. The American Astronomical Society. All rights reserved. Printed in the U.S.A. 


doi: 10.1 088/0004-637X/7 81/2/1 07 


FIRST HARD X-RAY DETECTION OF THE NON-THERMAL EMISSION AROUND THE ARCHES CLUSTER: 

MORPHOLOGY AND SPECTRAL STUDIES WITH NuSTAR 

Roman A. Kriyonos 1 , John A. Tomsick 1 , Franz E. Bauer 2,3 , Frederick K. Baganoff 4 , Nicolas M. Barriere 1 , 
Arash Bodaghee 1 , Steven E. Boggs 1 , Finn E. Christensen 5 * , William W. Craig 1,6 , Brian W. Grefenstette 7 , 
Charles J. Hailey 8 , Fiona A. Harrison 7 , JaeSub Hong 9 , Kristin K. Madsen 7 , Kaya Mori 8 , Melania Nynka 8 , 

Daniel Stern 10 , and William W. Zhang 11 
1 Space Science Lab, University of California, Berkeley, CA 94720, USA; krivonos@ssl.berkeley.edu 
2 Instituto de Astrofisica, Facultad de Flsica, Pontificia Universidad Catlica de Chile, Casilla 306, Santiago 22, Chile 
3 Space Science Institute, 4750 Walnut Street, Suite 205, Boulder, CO 80301, USA 
4 MIT Kavli Institute for Astrophysics and Space Research, Cambridge, MA 02139, USA 
5 DTU Space-National Space Institute, Technical University of Denmark, Elektrovej 327, DK-2800 Lyngby, Denmark 
6 Lawrence Livermore National Laboratory, Livermore, CA 94550, USA 
7 Cahill Center for Astronomy and Astrophysics, California Institute of Technology, Pasadena, CA 91125, USA 
8 Columbia Astrophysics Laboratory, Columbia University, New York, NY 10027, USA 
9 Harvard- Smithsonian Center for Astrophysics, 60 Garden Street, Cambridge, MA 02138, USA 
10 Jet Propulsion Laboratory, California Institute of Technology, Pasadena, CA 91109, USA 
11 NASA Goddard Space Flight Center, Greenbelt, MD 20771, USA 
Received 2013 September 30; accepted 2013 December 7; published 2014 January 15 


ABSTRACT 

The Arches cluster is a young, densely packed massive star cluster in our Galaxy that shows a high level of star 
formation activity. The nature of the extended non-thermal X-ray emission around the cluster remains unclear. The 
observed bright Fe Ka line emission at 6.4 keV from material that is neutral or in a low ionization state can be 
produced either by X-ray photoionization or by cosmic-ray particle bombardment or both. In this paper, we report 
on the first detection of the extended emission around the Arches cluster above 10 keV with the NuSTAR mission, 
and present results on its morphology and spectrum. The spatial distribution of the hard X-ray emission is found to 
be consistent with the broad region around the cluster where the 6.4 keV line is observed. The interpretation of the 
hard X-ray emission within the context of the X-ray reflection model puts a strong constraint on the luminosity of 
the possible illuminating hard X-ray source. The properties of the observed emission are also in broad agreement 
with the low-energy cosmic-ray proton excitation scenario. 
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1. INTRODUCTION 

The Arches cluster (GO. 12+0.02; Cotera et al. 1996; Serabyn 

et al. 1998) is a massive star cluster with a core that is about 

9" (~0.35 pc at 8 kpc) in radius (Figer et al. 1999), with an 

average mass density of ~3 x 10 5 M 0 pc -3 , containing more 
than 160 O-type stars (Figer et al. 2002). The Arches cluster 
is located in the inner Galactic center (GC) region at the 
projected angular distance of IT from the dynamic center of 
the Galaxy — Sagittarius A* (Sgr A*). 

The first serendipitous Chandra observations of the Arches 

cluster region revealed two bright sources (Al and A2) located 
in the cluster’s core which are surrounded by the cluster’s diffuse 
X-ray emission (Yusef-Zadeh et al. 2002). Al was later resolved 

into two distinct bright sources: AIN and A1S by Law & Yusef- 
Zadeh (2004) and Wang et al. (2006). Subsequent dedicated 

campaigns revealed a complicated picture of the Arches cluster 
X-ray emission, using spectral and morphological studies to 
establish the presence of thermal and non-thermal emission 
components. The thermal emission is thought to originate 
from multiple collisions between strong winds of massive 

stars (Chlebowski & Garmany 1991; Yusef-Zadeh et al. 2002; 
Wang et al. 2006; Capelli et al. 2011a). Diffuse non-thermal 
X-ray emission has been detected from a broad region around 
the cluster (Wang et al. 2006; Tsujimoto et al. 2007; Capelli 
et al. 2011b; Tatischeff et al. 2012). The non-thermal nature 

of this extended radiation is revealed by its bright fluorescent 


Fe Ka 6.4 keV line emission. The Fe Ka line at 6.4 keV, 
which is due to Fe that is neutral or in a low ionization 
state, is produced by ejecting a K-shell electron, either by hard 
X-ray photoionization from an external X-ray source, or by the 
collisional ionization induced by cosmic-ray (CR)- accelerated 
electrons or ions (mainly protons and a -particles). 

The mechanism producing the fluorescent line emission 
around the Arches cluster is not completely understood. The 
fluorescence may be the result of irradiation of the cloud by 
hard X-ray photons with energies above 7.1 keV (the K-edge 
of neutral iron). The required X-ray illumination might be 
associated with a nearby transient X-ray source (e.g., Churazov 
etal. 1993, for 1E1740.7-2942) or with past activity of Sgr A*, as 
suggested by Sunyaev et al. (1993) to explain the fluorescent line 
emission of the giant molecular cloud Sgr B2 in the GC region 
(Koyama et al. 1996; Sunyaev & Churazov 1998; Murakami 
et al. 2000; Revnivtsev et al. 2004a; Terrier et al. 2010). 
Alternatively, the excitation of cold material might be through 
the collisional ionization by CR particles (e.g., Predehl et al. 
2003; Yusef-Zadeh et al. 2007; Dogiel et al. 2009). Tatischeff 
et al. (2012, hereafter T12) argued that bright 6.4 keV line 
structures observed around the Arches cluster are very likely 
produced by bombardment of molecular gas by energetic ions, 
and they disfavored the possibility of accelerated electrons. T12 
suggest that the required large flux of low-energy cosmic ray 
(LECR) ions could be produced in the ongoing supersonic 
collision between the star cluster and a nearby molecular cloud. 
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Table 1 

Summary of the NuSTAR Observations Available for the Arches Cluster 


Observation 

Date 

Exposure a 

(ks) 

Off-axis b 

ID 

2012 

FPMA 

FPMB 

(arcmin) 

40010003001 

Oct 14 

24.2 


3.7 

40010005001 

Oct 15 

(17.0) 

26.0 

25.6 

4.1 

40010006001 

Oct 16 

(15.7) 

23.6 

(20.0) 

23.1 

7.5 



(12.0) 

(12.1) 



Notes. 

a Contains dead-time corrected nominal exposure, whereas the effective expo- 
sure, estimated for vignetting at 8 keV, is listed in parentheses. 
b The off-axis angle is calculated for the Arches cluster position with respect to 
the optical axis position on the detector. 

Previously, the lack of spectral and morphological studies 
above 10 keV precluded us from building a complete picture 
of the Arches cluster non-thermal emission and from making 
a definitive conclusion about the ionizing mechanism. With 
the first focused observations of the Arches cluster in the hard 
X-ray domain, we can shed light on the cause of the non-thermal 
X-ray emission detected around the cluster. 

We assume a distance of 8 kpc to the GC (Gillessen et al. 
2009) and the Arches cluster. Section 2 of the paper describes the 
Nuclear Spectroscopic Telescope Array ( NuSTAR ) observations 
and data reduction. In Section 3, we study the morphology of 
the Arches cluster using two-dimensional (2D) image fitting 
procedures and wavelet image decomposition techniques. In 
Section 4, we present a spectral analysis using the background 
approach specially developed for NuSTAR observations. The 
main results are discussed in Section 5 and summarized in 
Section 6. 

2. OBSERVATIONS AND DATA REDUCTION 

We observed Arches cluster during the GC region campaign 
with the NuSTAR (Harrison et al. 2013) in 2012 October. 
NuSTAR operates in wide energy band from 3 to 79 keV, 
carrying two identical co-aligned X-ray telescopes with an 
angular resolution of 18" (FWHM). The focal planes of each 
telescope, referred to as focal plane module A and B (FPMA 
and FPMB), provide a spectral resolution of 400 eV (FWHM) 
at 10 keV. 

The first part of the survey was carried out in the inner GC 
region covering the area between Sgr A* and the persistent 
low-mass X-ray binary IE 1743.1-2843 located 19' away 
(Porquet et al. 2003b; Del Santo et al. 2006). The Arches cluster 
serendipitously appears in the three NuSTAR observations listed 
in Table 1 . 

We processed the data from both modules using the NuSTAR 
Data Analysis Software (NuSTARDAS) vl.2.0, the 2013 June 
28 version of the NuSTAR Calibration Database (CAEDB), and 
HEASOFT v6.13. The data were filtered for periods of Earth 
occultation, high instrumental background due to South Atlantic 
Anomaly passages, and known bad/noisy detector pixels. 

The NuSTAR detectors are not entirely shielded from X-rays 
that pass outside the X-ray optics modules and fall directly 
onto the detectors. Such an unfocused flux or “stray-light” can 
be significant if there are bright X-ray sources within 2-3 
deg of the NuSTAR field of view. Stray-light from a point 
source has a uniform distribution on the detector and can be 


easily removed based on a geometrical model of the telescopes. 
We masked out the stray-light region of the FPMB detector 
(FPMA is clean) caused by two bright sources: GX 3+1 and 
the transient Swift J1745 10.8-262411 (see, e.g., Vovk et al. 
2012 and references herein), which was active at the time of 
the NuSTAR observations. The contaminated detector pixels 
have been flagged as “bad” and removed from the analysis 
at the stage of data screening in the NuSTARDAS pipeline run. 
The observation 40010003001 (FPMB) has been completely 
disregarded because stray-light completely covered the Arches 
cluster region. 

Stray-light photons from the extended sources and known 
celestial X-ray backgrounds cannot be so easily removed be- 
cause they fill the whole detector area and have a non-uniform 
pattern depending on the emission morphology. We treated this 
background differently for the imaging and spectral analysis. 
Regarding the former, the non-uniform background of the indi- 
vidual observations is effectively averaged when it is combined 
into the mosaic, and therefore, can be considered as nearly flat 
(however, see the discussion in Section 3.1). In contrast, the 
non-uniformity of the background has a noticeable effect on the 
spectral analysis since it is done for individual observations. 
To overcome this issue we applied the approach described in 
Section 4. 

After running the pipeline and screening the data, we end 
up with event lists for five NuSTAR observations containing the 
Arches cluster at different off-axis angles (Table 1). Note that 
the source appears at large off-axis distances, which causes the 
reduction of the efficiency (vignetting) and strong Point Spread 
Function (PSF) distortion. The total nominal (effective) dead- 
time corrected exposure time of 122.5 (76.7) ks comprises 73.8 
(44.7) ks and 48.7 (32.1) ks for FPMA and FPMB, respectively. 

3. MORPHOFOGY 

The data from both detector modules were used to construct 
the mosaic images (Figure 1) around the Arches cluster in 
different energy bands. We used the 3-10 keV NuSTAR band 
as a proxy for the standard 2-10 keV energy range; 6-7 keV 
band as containing Fe Ka emission lines at 6.4 and 6.7 keV, 
which are unresolved by NuSTAR ; and 10-20 keV as a hard 
X-ray band specific to NuSTAR. The mosaics in the energy bands 
3-79, 20-40, and 40-79 keV are also shown in Figure 1. For 
each image an average (local) background level was determined 
within an annular region centered on the Arches cluster, as 
defined in Table 2 and illustrated in Figure 2. The background 
estimate was subsequently subtracted from each individual 
image, prior to the construction of the image mosaic. This is 
done because FPMB has a higher background count rate than 
FPMA, requiring that we adjust the relative detector background 
normalizations between modules. Note that the visible gradient 
in Figure 2 is produced mainly by the emission from the GC 
molecular clouds (see, e.g., Ponti et al. 2010). 

The exposure map, utilized for mosaic creation, has not been 
corrected for vignetting effects. We checked that the vignetting 
correction does not significantly improve the results of the 
imaging analysis. Hereafter, we assume no vignetting correction 
unless otherwise stated. 

The primary goal of the current paper is to study the non- 
thermal emission surrounding the Arches cluster. The emission 
is revealed by fluorescent Fe Ka 6.4 keV line, as has been 
observed by XMM-Newton and reported in the recent papers 
by Capelli et al. (2011a) and T12. We adopted T12 sky 
regions describing the core of the Arches cluster prominent 
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Figure 1. NuSTAR mosaic images of the Arches cluster in different energy bands. The images are background-subtracted, and adaptively smoothed with task 
dmimgadapt from Cl AO-4. 5 using a tophat kernel. The region used for local background estimation is shown in Figure 2. The cyan circle (in black on the two upper 
images) denotes the Arches cluster emission prominent at 6.7 keV (a He-like Fe Ka line), and this region is denoted as “Cluster” in Table 2. The green ellipse represents 
the surrounding region characterized by the bright 6.4 keV line emission (“Cloud” region in Table 2). Both regions have been adopted from T12 for reference. The 
images are shown in the same log-scale color map shown on the right. The color bar is in units of counts s -1 pixel -1 , and the angular size of the image pixel is 2". 
The white grid indicates Galactic coordinates in degrees. The compass sign shows the image alignment in equatorial coordinates, with north up and east to the left. 
The image orientation remains unchanged throughout the paper. 

(A color version of this figure is available in the online journal.) 


Table 2 

Definitions of the Sky Regions 


R.A. (J2000) 

Decl. (J2000) 

Parameters 

17 h 45 m 5K3 

Cluster (circle) 
—28° 49' 19" 

15" 

17 h 45 m 51)0 

Cloud (ellipse) 
-28° 49' 16" 

25", 59", 155° 

17 h 45 m 5K3 

—28° 49' 19" 

15" (excl.) 

17 h 45 m 50?3 

Background (annulus) 
—28° 49' 19" 

185" 

17 h 45 m 5(K3 

-28° 49' 19" 

100" (excl.) 


Notes. Central position and radius for circular regions, and semi- 
minor/major axes and rotation angle for the elliptical regions. The 
rotation is defined counter clockwise relative to North (upward). 
“Cluster” and “Cloud” regions are illustrated in Figure 1. The 
excluded regions are marked as “excl.” 


in Fe Ka 6.7 keV line (“Cluster” region) and the surrounding 
region characterized by the fluorescent 6.4 keV line (called 
“Cloud” region). The 6.7 keV line, attributed to the cluster core, 
originates from the hot, thermally ionized plasma, probably 
from the colliding winds of the massive stars. The 6.4 keV line 
emission clearly extends beyond the Arches cluster core, and 
comes from material in neutral or low ionization states filling 
(in projection) an extended area represented by an ellipse as 
illustrated in Figure 1. The ellipse is oriented approximately 
in the north-south direction with increased brightness to the 
south. The parameters of the background and reference regions 
are listed in Table 2. 



Figure 2. Definition of the local background region around the Arches cluster. 
The image shows the NuSTAR mosaic of the Arches cluster in the 3-10 keV 
energy band (see Figure 1 for the reference). The yellow annulus denotes the 
local background region used in the mosaic build up and fitting procedure 
(Section 3.1). The parameters of the annulus are listed in Table 2. 

(A color version of this figure is available in the online journal.) 


3.1. Cluster Emission 

The Arches cluster has been relatively well studied in the 
standard 2-10 keV energy band with the Chandra and XMM- 
Newton observatories. Close attention has been paid to the 
central part of the cluster, where thermal emission, thought to be 
from colliding winds of giant stars, is observed. With data from 
NuSTAR , we can obtain new insight into the hard X-ray emission 
of the Arches cluster. However, even with angular resolution that 
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Figure 3. Average NuSTAR PSF in the 3-10 keV energy range. The blue 
contours denote PSF encircled energy levels for 10%, 30%, 50%, 80%, and 
90%, respectively, going from the center outwards. The PSF was averaged over 
five observations listed in Table 1 . The green ellipse and black circle represent 
Arches cluster emission regions; see Figure 1 for details. 

(A color version of this figure is available in the online journal.) 


is unprecedented for hard X-ray astronomy (18" FWHM), it is 
challenging to resolve emission in the core, which has a radius 
of -15". 

The 3-10 keV mosaic (Figure 1) demonstrates a strong excess 
in the core of the Arches cluster and apparent extended emission 
around it. To decouple these emission components and find the 
centroid position of the core, we fitted a spatial model convolved 
with an averaged PSF on the combined data set using the Sherpa 
package (Freeman et al. 2001), which is part of the CIAO-4.5 12 
software (Fruscione et al. 2006). 

To make a PSF function representing the average over the 
different source offsets in the ObsIDs, we used the inflight- 
calibrated PSF shape stored in the CALDB. The PSF images 
were weighted according to the exposure time of each data set. 
The resulting averaged PSF is shown in Figure 3. Its shape is 
strongly distorted because of the large off-axis angles (on-axis 
PSF would be symmetric). Importantly, the PSF elongation is 
not co-aligned with the 6.4 keV line cloud region. For the fitting 
procedure, we assume no significant change of PSF shape inside 
the radius of 60" around the Arches cluster (see below). Note 
that 90% of the averaged PSF encircled energy is situated at 
or inside of the inner radius of the annular region used for the 
background estimation, depending on the adopted energy used. 
Using this fact, we estimated the background count rate in the 
separate fitting procedure without PSF convolution. 

Given the spatial extent of the emission around the Arches 
Cluster as found by Chandra and XMM-Newton , we want 
to confirm whether we can resolve this same emission with 
NuSTAR. This will give us confidence that we may also be able 
to spatially resolve the Arches cluster above 10 keV. 

We represented the Arches cluster core emission by a 2D 
Gaussian. Due to the low count rate per pixel, we used Poisson 
statistics in the fitting procedure. The background term has been 
constrained in the annulus region shown in Figure 2. After 
running the fitting procedure with the 2D Gaussian model, we 
found that the 3-10 keV emission from the Arches cluster was 


12 http://cxc.cfa.harvard.edu/ciao 



Figure 4. Radial profile of the Arches cluster region in the 3-10 keV energy 
band. Lines show the best-fit model and its components convolved with the 
averaged PSF. 

(A color version of this figure is available in the online journal.) 


Table 3 

Best-fit Model Parameters for the PSF Fitting Procedure of the 3-10 keV 
Arches Cluster Core and Halo Emissions 


Parameter 

Value 

Norm, (x 10 -5 ) 

Background (constant) 

1.18 ± 0.01 counts s -1 pixel -1 

Center R.A. a 

Core (2D Gaussian) 

17 h 50 m 5K43 

Center Decl. a 

— 28°49'23"07 

FWHM 

470 (fixed) 

Norm, (x 10 -3 ) 

2.2 ± 0.3 counts s -1 pixel -1 

Center R.A. a 

Halo (2D Gaussian) 

17 h 45 m 5K62 

Center Decl. a 

-28°49'47717 

FWHM 

72'.'4 +6 >7 

Norm, (x 10 -5 ) 

-4.7 

1 -77+‘o 4 37 counts s -1 pixel -1 

Cash statistic 

4661.91 

Degrees of freedom 

2820 


Note. a J2000 epoch, the Icr positional uncertainty is less than 2" and 6"-8" 
for the core and halo components, respectively; see 2D confidence regions in 
Figure 6 for details. Model normalization factor is marked as “Norm.” 


not acceptably fitted with a single point source, and an extended 
“halo” around the cluster had to be considered. This finding 
is in line with the detection of extended emission surrounding 
the Arches cluster as observed by Chandra and XMM-Newton 
(Yusef-Zadeh et al. 2002; Law & Yusef-Zadeh 2004; Wang et al. 
2006). 

We constrained the fitting area to a circle centered on the 
core and extending to a radius of 60". The model includes (1) 
a narrow 2D Gaussian of fixed width (4" FWHM, chosen to 
emulate the PSF smearing effect caused by incomplete aspect 
reconstruction), representing the cluster’s core, (2) a wide 2D 
Gaussian describing the halo, and (3) a constant background 
term. Note that the visible in Figure 2 background gradient has 
no significant effect on the fitting procedure. The best-fit model 
parameters are listed in Table 3. Figure 4 shows the radial profile 
of the Arches cluster intensity distribution centered at the best-fit 
position of the core. 

To demonstrate a surface brightness distribution of the Arches 
cluster halo, we show in Figure 5 the residuals after subtracting 
off a best-fit point source corresponding to the core emission. 
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Figure 5. NuSTAR 3-10 keV residual image of the Arches cluster after removal 
of the emission from the cluster’s core and background subtraction. The details 
of the region description can be found in Figure 1. The color bar units are 
counts s -1 pixel -1 and the angular pixel size is 2". 

(A color version of this figure is available in the online journal.) 

The position of the cluster’s core is spatially consistent with 
the brightest source, Al, detected by Chandra in the Arches 
cluster (Yusef-Zadeh et al. 2002) and later resolved into the two 
sources A1S/N (Law & Yusef-Zadeh 2004; Wang et al. 2006). 
Figure 6 shows confidence regions for the position uncertainty 
of the core and halo components. 

Note that due to the NuSTAR optics design, the celestial 
coordinates of each incoming photon are determined during 
the post-processing. This task is complicated by distortions 
due to thermal bending and external forces acting on the 
mast during orbit. The systematic offset of the source position 
determination can be as high as 10". However, the fact that 
the position of the Arches cluster core emission detected with 
NuSTAR is in agreement (within the angular separation of ~ 1"5) 
with the catalogued coordinates of the two bright Chandra 
sources A1S/N, tells us that the described systematic offsets 
are negligible for these observations. 


Our ability to centroid a weak source with NuSTAR is also 
subject to uncertainties based on how positions are probabilis- 
tically assigned to photons in the NuSTARDAS pipeline and 
the finite spatial resolution of the NuSTAR pixels. To estimate 
any scatter introduced by this process, we ran the NuSTARDAS 
pipeline 50 times, repeated out analyses, and found that the sky 
position of the Gaussian core has an rms scatter of O'.' 9. We adopt 
this as the systematic error on the cluster core position. 

The width of the 2D Gaussian halo model component does 
not have strong physical meaning since this component is used 
specifically to describe local background of the core. However, 
we should note that the center position of the halo component 
does not coincide with the core centroid, and the trend suggests 
that the halo centroid falls outside the Arches cluster (Figure 6). 


3.2. Wavelet Decomposition 

We studied in Section 3.1 the morphology of the Arches 
cluster in the 3-10 keV band fitting data with a parametric 
model. We constrained the central part of the cluster, determined 
its centroid position, flux, and revealed strong excess outside 
the cluster core. Since we observe extended emission with a 
complicated projected morphology, it becomes quite difficult 
to represent this emission with an analytical model, i.e., to 
parameterize it. In this case, non-par ametric methods can be 
used for morphological studies. One such method is a wavelet 
image transform. 

For the decomposition of NuSTAR images in different energy 
bands, we used the so called a trous digital wavelet transform 
algorithm (Starck & Murtagh 1994; Slezak et al. 1994; Vikhlinin 
et al. 1997; Krivonos et al. 2010). The main advantage of 
this algorithm over other wavelet methods is the fact that it 
decomposes the image into a linear combination of layers 
or “scales,” each of them emphasizes the structures with a 
characteristic size of ~2 /_1 pixels, which is 2 J arcsec in the 
case of our NuSTAR mosaic image resolution (1 pixel = 2"). 
Thus, low (high) J values correspond to small (large) spatial 
variations. We used the WVDECOMP wavelet decomposition 
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Figure 6. Left: Arches cluster core and halo centroid determination on the 3-10 keV NuSTAR mosaic image (use Figure 1 for the reference). Blue contours show 
lor, 2 or, and 3cr confidence intervals for the core (solid blue) and the halo (dashed blue) centroids. The centroid positions of the core and halo emission are marked 
as blue crosses. Right: a close-up of the region of the Arches core. Black points mark positions of three bright sources A1S/N and A2 that are resolved by Chandra 
observations of the Arches cluster (Wang et al. 2006). Note that a systematic scatter of O'/ 9 (rms) is not taken into account (see Section 3.1). The color bar units are 
counts s -1 pixel -1 and the pixel size is 2". 

(A color version of this figure is available in the online journal.) 
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Figure 7. Wavelet decomposition of NuSTAR mosaic images of the Arches cluster in 3-10, 6-7, 10-20, and 20-40 keV energy bands. All structures with detection 
significance below 5 a have been removed from the images. The units of the color bar on the right are counts s -1 pixel -1 , where the pixel size is 2" . 

(A color version of this figure is available in the online journal.) 



Figure 8. Stacked wavelet layers (see Figure 7) of the Arches cluster NuSTAR mosaic images in the (from left to right) 3-10, 6-7, and 10-20 keV energy bands. All 
structures with detection significance below 5cr have been removed from the images. The units of the images are counts s -1 pixel -1 and the pixel size is 2". 

(A color version of this figure is available in the online journal.) 


code developed by Vikhlinin et al. (1995, 1998) as a part of 
ZHTOOLS software package. 13 

We ran WVDECOMP on mosaic images for scales 1-7 
(2"- 128"), detecting all significant structures above the 5a level. 
To keep the photon statistics on the images, we provided counts 
and exposure maps to WVDECOMP. As mentioned in Section 3, 
FPMB has a higher background count rate than FPMA. In fact, 
this can mimic extended structures on the exposure-corrected 
image. To diminish this effect, we built the background map 
for a given observation by multiplying the exposure map and 
background count rate measured in the annulus region shown in 
Figure 2. The individual background maps were combined into 
the mosaic and used as a part of the WVDECOMP algorithm. 

Figure 7 shows wavelet scales 3-6 of the Arches cluster 
mosaic images in 3-10, 6-7, and 10-20 keV. The lowest scales 


1 3 http : //hea- www. harvard . edu/RD/ zhtools 


(1,2) do not contain significantly detected (5a) structures. 
The highest scale (7) contains large-scale diffuse background 
attributed to the Galactic center molecular complex zone (Ponti 
et al. 2010) and Galactic Ridge X-ray Emission (GRXE; see, 
e.g., Revnivtsev et al. 2006). We combined scales 3-6, which 
only contain structures related to the Arches cluster into the 
mosaics shown in Figure 8. The main results from the WT 
decomposition can be summarized as follows (refer to Figures 7 
and 8). 

3-10 keV. The core of the Arches cluster is clearly visible 
on the scales 7 = 3-4 and the halo emission is revealed on 
scale 7 = 5. The apparent position of the halo is consistent with 
fitting results in Section 3.1 showing a shifted position with 
respect to the cluster core. 

6-7 keV. This energy band contains two blended Fe Ka lines 
at 6.4 and 6.7 keV. We speculate that these lines appear at scales 
7=4 (6.7 keV) and 7 = 5 (6.4 keV). The 6.7 keV line is 
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spatially consistent with core emission, while the 6.4 keV line 
emission follows the cloud region as expected. These results 
are in full agreement with XMM-Newton observations (T12). 
In contrast to the 3-10 keV band, the 6.7 keV line emission 
distribution does not demonstrate high-frequency detection at 
7 = 3, which may indicate its broader spatial distribution 
in comparison with the 3-10 keV continuum emission of the 
Arches cluster. However, we should note that the 6-7 keV 
mosaic image is not continuum subtracted. 

10-20 keV. The Arches cluster does not show any signifi- 
cant (>5cr) point-like structures in the core region, which is 
expected due to the rapid fall of the ~1.7 keV thermal emis- 
sion above 10 keV. On the other hand, the hard X-ray emis- 
sion is significantly detected at scales 7 = 5-6. Note that 
the observed extended structure cannot be caused by the wide 
wings of the PSF — as is demonstrated in Figure 3, the PSF 
shape is not aligned with the cloud region morphology. The 
observed 10-20 keV continuum emission apparently coincides 
with Fe Ka fluorescent 6.4 keV line emission observed with 
XMM-Newton (T12). According to T12, the 6.4 keV line emis- 
sion demonstrates brightening in the south-east part of the 
ellipse, which is also noticeable at scale 7 = 5. 

20-40 keV. The mosaic map of the Arches cluster shown in 
Figure 1 does not allow to make a definitive conclusion about 
the source detection. However, wavelet decomposition shows 
weak excess at scale 7 = 6. This fact is implicitly confirmed by 
detection of the emission up to 40 keV in the spectral analysis of 
Section 4. The apparent hard X-ray emission is located outside 
the cluster and consistent with general extent of 10-20 keV 
emission. 

4. SPECTROSCOPY 

We extracted spectral information for the Arches cluster from 
the circular region with radius of 50" centered at the clus- 
ter’s core (Table 2), including the core and extended emis- 
sion. The spectra have been extracted from all observations 
using the nuproducts task applying PSF and ghost ray correc- 
tions to the effective area Ancillary Response File. 

The extraction of the corresponding background spectra is 
complicated by the strong CXB and GRXE stray-light illumi- 
nation causing a non-uniform pattern on the detectors with a 
gradient of a factor of two. For that reason, the background 
spectrum extracted from the source-free region over- or under- 
predicts the background count rate at the position of the source. 
This difficulty (not relevant for spectra of bright sources) causes 
strong systematic noise for the background- subtracted spectrum 
of weak sources such as Arches. To overcome this, we use the 
fact that the stray-light background does not significantly change 
for sky pointings separated by only l(y-2(y, and therefore we 
utilize nearby observations. Note that this method differs from 
that used in the imaging analysis. In the paired observation, we 
extract the background spectrum from the same detector region 
where the Arches cluster was observed. Obviously, this region 
should not contain any other diffuse or point-like source or stray- 
light from nearby bright sources. Table 4 lists the observations 
used for source and background spectrum extraction. For ObsID 
40010006001, we had three observations which were suitable 
for background spectrum extractions and thus averaged together 
for better statistics. 

Using the above-mentioned method, we have constructed 
one FPMB and three FPMA spectra. The latter have been 
further grouped together using the FTOOL task MATHPHA. 
The corresponding background spectra and responses have 


Table 4 

Observations Used for Source and Background Spectral Extraction 


Source Spectrum 
(ObsID a ) 

Background Spectrum 
(ObsID a ) 

4001 000300 lA b 

40010001001A 

400 10005001 A 

4001 000300 lA b 

4001 000600 1A/B C 

4001 000 1 002 A/B 


4001 000300 1A/B 


400 10005001 A/B 


Notes. 

a The name of the focal plane module (A or B) is attached to 
the end of the ObsID. 

b The corresponding FPMB data has strong stray-light 
contamination, see Section 2. 

c The background spectrum has been averaged over three 
observations. 


been averaged accordingly. Finally, we have grouped spectrum 
channels to have at least 30 counts bin -1 using GRPPHA tool. 
Figure 9 shows the FPMA and FPMB spectra of the Arches 
cluster. 

According to the morphological studies presented in Sec- 
tion 3, we attribute the hard spectral component above 10 keV 
to the fluorescing cloud around the Arches cluster. 

We first approximated the observed spectra with the colli- 
sionally ionized plasma emission model APEC representing the 
cluster’s core thermal emission, and a non-thermal power-law 
continuum with a Gaussian line at 6.4 keV, modeling the cloud 
emission outside the cluster (Figure 9). Due to the fact that 
NuSTAR has a limited spectral resolution, the centroid energy 
of the Gaussian line and its width have been fixed at 6.4 keV 
and 0.1 keV, respectively. All the emission components were 
subject to a line-of-sight photoelectric absorption model WABS 
in XSPEC. The fitting procedure does not allow us to constrain 
the absorption, and we fixed it at Nu = 9.5 x 10 22 cm -2 mea- 
sured by T12 (also using WABS) in the core of the Arches 
cluster. The metallicity of the X-ray emitting plasma was not 
constrained either, and we fixed it to Z = 1 .7 Z© as per T 12. The 
best-fit results obtained with this model (referred to as model 
1 in the discussion) are shown in Table 5. All fitted parameters 
in this table and in the following discussion contain uncertainty 
estimation at the 90% confidence level. 

Model 1 provides a reliable fit to the data by constraining the 
temperature of the core emission, the power-law component and 
the intensity of the 6.4 keV line. The 1.76 +q 2 6 9 keV temperature 
of the thermal component is in agreement with that recently 
measured by T12 in the cluster’s core (1.79 +qq 6 5 keV from 
“model 1” in T12). The power-law slope T = 1.62 ± 0.31 
is consistent with T = 1.6+q 2 from T12, measured in the cloud 
region. T12 finds the slope and 6.4 keV line equivalent width 
(EW 6 . 4 kev) of the cloud region compatible with the LECR ion 
model, where collisional ionization of the cloud by CR protons 
produces the 6.4 keV fluorescence line. According to the XMM- 
Newton study of the Arches cluster by T12, the LECR electron 
model (LECRe) requires an ambient Fe abundance >3 times the 
solar value to account for the measured EW of 1.2 ± 0.2 keV, 
which makes the measured properties of the non-thermal cloud 
region emission hardly compatible with the LECR electron 
model. We consider only the LECR proton (LECR/?) model 
in the following analysis. 

We fitted the spectra with the XSPEC model WABS x 
(APEC + LECR/?), henceforth “model 2.” The model depends 


7 


The Astrophysical Journal, 781:107 (llpp), 2014 February 1 


Krivonos et al. 



1*1 

-iM MilM iiiiiliii!itfr^ 

FPMA 

FPMB 

_ 


| ■ -f- -j- 


5 10 20 50 

Energy (keV) 


> CO 
CD 1 “ 

tltAMkuJ 

FPMA 

FPMB 

05 7 “ 

'Tn 



10 

05 1 

O 

co 



CM - 

J- 

CO 



~5 10 20 50 ~ 

Energy (keV) 


Figure 9. X-ray spectrum of the Arches cluster as measured with NuSTAR FPMA (black) and FPMB (red). The spectra are fitted with models all containing APEC 
component (magenta) and: model 1 (left) — power-law (blue) and Gaussian (green); model 2 (middle) — LECRp (blue); model 3 (right) — REFLIONX (blue); see 
Table 5 for details. Solid lines represent model components for FPMA, and dotted for FPMB, respectively. 

(A color version of this figure is available in the online journal.) 


Table 5 


Spectral Analysis of the Arches Cluster X-Ray Emission. 

Parameter 

(Unit) 

Model 1 

Model 2 

Model 3 

C 


1.06 ±0.14 

1.05 ±0.14 

1.07 ±0.14 

N H 

(10 22 cm -2 ) 

9.5 (fixed) 

9.5 (fixed) 

9.5 (fixed) 

Z/Z 0 


1.7 (fixed) 

1.7 (fixed) 

1.7 (fixed) 

kT 

(keV) 

1 76 +0 - 36 
10 -0.29 

1 AQ+0.31 
A ‘ oy -0.26 

1 A4+0.37 
1,0 -0.23 

hr 

(see notes) 

36 -°-12.9 

39 6+ 16 - 7 
11.9 

43 2 +16 - 9 

Eg A keV 

(keV) 

6.4 (fixed) 



£6.4 keV 

(10 -5 ph cm -2 s -1 ) 

1.24 ±0.5 



r 


1.62 ±0.31 



fp.i. 

(10 -12 erg cm -2 s -1 over 3-20 keV) 

1 4Q+0.26 
A -^-0.24 



EW6.4keV 

(keV) 

1 l +0 ' 7 



A 

(H-atoms cm -2 ) 


5 x 10 24 (fixed) 


s 



1 A5+0.59 
A * o;5 -0.55 


£min 

(keV/n) 


10 4 (fixed) 


Alecr 

(10 -8 erg cm -2 s -1 ) 


8-08-1.62 


Erx 




2.93 ±0.31 

$ 

(erg cm s _1 ) 



10 (fixed) 

/rx 

( 10 - 5 ) 



6.09 ± 5.00 

Reduced x 2 (d.o.f.) 


0.91 (134) 

0.89 (136) 

0.91 (140) 


Notes. Model 1: C x WABS x (APEC + Gaussian + power law); Model 2: C x WABS x (APEC + LECR/?); Model 3: C x WABS x 
(APEC + REFLIONX). The cross-calibration constant term C is fixed at unity for FPMA data and fitted for FPMB. The absorption 
column density is noted as Ar. APEC thermal plasma model is characterized by kT, Z/Z Q , and Et parameters, respectively, describing 
temperature, metallicity relative to solar, and normalization in units of 1 0“ 1 8 f n e nndV /(4 tt D 2 ), where n e and izr are the electron and 
proton number densities in units of cm -3 , and D is the distance in cm. £ 5 . 4 keV and £6.4 keV- centroid energy and flux of the neutral or 
low-ionization Fe Ka line. T and 7 p j.: index and 3-20 keV pegged normalization of the power-law component (model PEGPWRLW). 
EW 6 .4kev: EW of the 6.4 keV line with respect to the power-law continuum. A, 5 , £ m in, and Alecr- LECR ions path length, source 
spectrum index, minimum energy, and model normalization. By definition d W/dt = 4nD 2 Alecr is the power injected in the interaction 
region by primary CR protons with energies between E m \ n and £ max = 1 GeV, where D is the distance to the source. Trx, and 7 rx: 
the power-law slope of the incident radiation, ionization parameter (defined as £ = 47r£ to t/^H where £ to t is the total illuminating flux, 
and «h is the density of the reflector), and normalization of REFLIONX model (Ross & Fabian 2005). The metallicity parameter Z for 
LECR and REFLIONX models is fixed at 1.7 Z Q . Reduced x 2 (d.o.f.): reduced x 2 and degrees of freedom. 


on the LECR path length in the ambient medium, A, the 
minimum energy of the CR particles entering the cloud, E^ n , 
the power-law index of the CR source energy spectrum, s 9 and 
the metallicity of the X-ray emission region, Z. We fixed A 
and E m { n parameters of the LECRp model according to T12, 
making free slope a and normalization Alecr- The best-fit results 
obtained with model 2 are given in Table 5 and the spectral 
components are shown in Figure 9. 

The A/wSTAA-measured spectral index of the injected low- 
energy protons s = 1.65+q 5 5 9 5 is in general agreement with 
the T12 value (1.9 + q 5 6 ). The best-fit LECR model normal- 


ization for the fixed E m \ n = 10 MeV nucleon 1 Alecr = 
(8.08+ 1 1 ,8 6 2 2 ) x 10 -8 erg cm -2 s -1 is somewhat higher than that 
measured by T12 (Alecr = (5.6+q 7 3 ) x 10 -8 erg cm -2 s -1 ). 
The corresponding power injected by LECR protons lies in the 
range (5-8) x 10 38 erg s -1 . 

To test the X-ray reflection scenario of the Arches cluster 
non-thermal emission, we applied the self-consistent reflec- 
tion model REFLIONX 14 in XSPEC describing the reflected 


14 http://heasarc.gsfc.nasa.gov/xanadu/xspec/models/reflion.html 
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spectrum for an optically thick atmosphere of constant density, 
illuminated by radiation with a power-law spectrum (Ross & 
Fabian 2005), henceforth “model 3.” The REFLIONX spectrum 
contains modeled fluorescence lines and continuum emission 
(Figure 9). Developed for the surface of hot accretion disks, 
REFLIONX can still be applied for cold material around the 
Arches cluster, fixing the ionization parameter § at the lowest 
allowed value of 10 erg cm s -1 . The best-fit parameters of model 
3 are shown in Table 5. 

5. DISCUSSION 

In this section, we discuss the results of this work in the 
context of two scenarios proposed to explain the excitation 
mechanism of the fluorescent emission around the Arches 
cluster. 

5.1. CR Bombardment Hypothesis 

A particle bombardment excitation scenario for the Fe fluo- 
rescence line emission involves the interaction of CR ions with 
cold molecular gas. T12 developed a steady- state slab model in 
which accelerated particles penetrate a cloud of neutral gas at a 
constant rate until they either stop or escape from the cloud. The 
model describes the properties of the neutral Fe Ka line emis- 
sion excited by impacts of LECR ions. These authors argue that 
fluorescent 6.4 keV line emission of the Arches cluster is ex- 
cited by LECR protons with a relatively soft incident spectrum 
s = 1.9 !°o 5 6 . Our NuSTAR spectral analysis (model 2) implies 
a somewhat harder slope a = 1.65t° 0 5 5 9 5 , however, both mea- 
surements are consistent within the error intervals. T12 notes 
that production of 6.4 keV line photons by LECR collisional 
ionization is relatively inefficient and that high power in LECRs 
is needed to produce the observed fluorescent Fe Ka line flux. 
As the authors suggest, a large flux of LECR ions could be 
produced in the ongoing supersonic collision between the star 
cluster and nearby molecular cloud. 

Note that LECR/? model provides an equally good fit as the 
reflection model discussed below. 

5.2. X-Ray Reflection Hypothesis 

By definition, the excitation of the neutral material of the 
cloud by hard X-ray irradiation implies a transient or persistent 
X-ray source outside or within the cloud. 

Assuming that the irradiation source is located in the Arches 
cluster, the source luminosity required to produce the observed 
6.4 keV line flux can be estimated from the fluorescent model de- 
scribed in Sunyaev & Churazov (1998). Applying this approach 
on the XMM-Newton data, T12 showed that the luminosity of 
the source in the cluster must be at least 10 36 erg s _1 . However, 
the persistent X-ray luminosity is far below this level. Capelli 
et al. (2011a) recently reported X-ray flaring activity in the core 
of the Arches cluster, probably originating from stellar wind in- 
teractions in one or more massive binary systems. Such transient 
activity, in principle, can trigger fluorescent emission of the neu- 
tral material outside the Arches cluster, but the peak luminosity 
of the observed flare is far below the threshold mentioned above. 

The illuminating source could also be a transient X-ray binary 
system. Their typical X-ray luminosity can reach levels on the 
order of 10 36 -10 38 erg s -1 for high- and low-mass systems. 
If the bright state of an X-ray binary lasted for more than 
8 yr at or above 10 37 erg s -1 , it would be sufficient for the 
observed 6.4 keV flux (Capelli et al. 2011b). However, such a 


long outburst is much longer than typical flaring timescales of 
X-ray transients (Degenaar & Wijnands 2010), though, at least 
one exceptional case is known, the neutron star X-ray transient 
KS 1731-260 (Sunyaev & Kwant Team 1989; Sunyaev et al. 
1990) which became quiescent after more than a decade in the 
bright state (Wijnands et al. 2001 ; Revnivtsev & Sunyaev 2002). 

The AwSTAR-measured EW of the 6.4 keV line with respect 
to the power-law continuum, EW 6 . 4 k e v = LU° 0 7 5 keV, is 
in agreement with previous studies performed on spatially 
resolved fluorescent clouds around the Arches cluster (Capelli 
et al. 2011b, T12). In the context of an X-ray reflection 
scenario, the relatively high EW of the Fe Ka line indicates a 
scattered emission component (assuming normal Fe abundance) 
or strong photoabsorption (Sunyaev & Churazov 1998), which 
can explain the fact that the direct source of photoionization is 
not seen by the observer. However, according to NuSTAR maps 
of the Arches cluster above 10 keV (Figure 1), for which the 
efficiency of photoabsorption rapidly falls, no significant source 
has been detected. 

The Arches cluster spectrum is well-fitted with a model 
composed of an optically thin thermal plasma and an 
X-ray reflection component (Figure 9). The spectrum of the illu- 
minating X-ray source is estimated to have relatively soft slope 
r RX = 2.93 ±0.31. We may assess the luminosity of the primary 
source L RX required by the reflection component following the 
method used by Reis et al. (2013) to estimate the reflection frac- 
tion R. This parameter is defined as a measure of the ratio of the 
reflected flux to the continuum power-law flux. REFLIONX 
calculates R by dividing the extrapolated (1 eV-1000 keV) 
REFLIONX flux by the 100 eV-1000 keV power-law flux. On 
the other hand, the reflection scaling factor R is in fact the solid 
angle R = £1/2tx subtended by the reflecting cloud (R = 0 
means no reflection). We estimated the unabsorbed 0.3-8 keV 
luminosity L RX (at 8 kpc) for the tabulated R values — 0.5, 1.0 
(source above a plane), 1.5, and 2.0 (source inside the cloud) as 
2.5, 1.3, 0.8 and 0.6 x 10 36 erg s _1 , respectively. Wang et al. 
(2006) reported the 0.3-8 keV luminosity of the bright Arches 
cluster sources AIN/S and A2, respectively, as 7.2 (AIN), 11 
(A1S), and 4.6 x 10 33 erg s -1 , which is 1-2 orders of magnitude 
below the luminosity of the illuminating continuum required by 
the reflected component. This constraint is consistent with what 
can be drawn from the observed 6.4 keV line flux, using the 
general idea of Sunyaev & Churazov (1998), as demonstrated 
by T12. 

The interpretation of the observed spectra with the X-ray 
reflection model confirms that the Arches cluster itself is 
unlikely to be the source of the cloud fluorescence emission. 

The non-detection of Fe Ka line flux variability over more 
than 8 yr of XMM-Newton observations (Capelli et al. 2011b, 
T12) puts strong limits on the photoionization excitation of 
molecular clouds around the Arches cluster by a hard X-ray 
source. Nevertheless, the plausibility of the scenario with an 
earlier outburst of Sgr A* is under debate. The large-scale spatial 
distribution and time variation of the Fe-K fluorescence seen in 
the GC region reveal the apparent superluminal propagation 
of a light front illuminating a molecular nebula (Ponti et al. 
2010). In this context, it is quite reasonable to anticipate 
X-ray photoionization of the molecular clouds near the Arches 
cluster by Sgr A*. A possible decrease in the 6.4 keV flux with 
time would be expected in this scenario. For instance, strong 
variability of the 6.4 keV emission has been detected in Sgr B2 
(Terrier et al. 2010) after 10 yr of a persistent state as shown 
by Revnivtsev et al. (2004a); however, see Inui et al. (2009). 
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Capelli et al. (201 lb) found the relatively constant 6.4 keV flux 
of the Arches cluster clouds consistent with a long outburst from 
Sgr A* at a constant luminosity of 1 .4 x 10 39 erg s -1 (Ponti et al. 
2010). On the other hand, Clavel et al. (2013) possibly revealed 
a highly variable active phase of Sgr A* during the past few 
hundred years. 

The measured power-law slope of the illuminating source 
in the model 3 (Section 4) is softer (T ~ 3) than observed 
for Sgr A* flares (Porquet et al. 2003a, 2008; Nowak et al. 
2012). Additionally, the slope is softer than that estimated 
by Revnivtsev et al. (2004b) for Sgr B2 in the context of 
Sgr A* reflected emission. However, we should note that 
REFLIONX calculates the reflected spectrum for an optically 
thick atmosphere while the molecular clouds close to the Arches 
cluster exhibit transparency at t < 1 (Capelli et al. 2011b). A 
more realistic reflection model at lower r values and different 
source-cloud geometry settings can produce a harder spectrum 
of the primary emission. 

6. SUMMARY 

In this paper, we present results on the first focused hard 
X-ray observation of the Arches cluster performed by NuSTAR. 
We report on significant detection of the hard X-ray emission 
around the cluster with a spectrum that extends for the first time 
up to ~40 keV. The main results can be summarized as follows. 

1. The NuSTAR mosaic image of the Arches cluster shows 
a core component and extended halo emission in the 
3-10 keV energy band. The centroid position of the core 
emission coincides with the bright sources AIN/S resolved 
with Chandra in the Arches cluster. The halo emission 
around the cluster also follows the morphology observed 
with Chandra. The centroid position of the halo emission 
does not coincide with the core and tends to be outside the 
Arches cluster. 

2. Wavelet decomposition of sky mosaics does not show any 
significant (5a) detection of point-like sources in the core of 
the Arches cluster above 10 keV. The continuum 10-20 keV 
emission is significantly detected well beyond the Arches 
cluster with a spatial morphology consistent with the 
Fe Ka fluorescent 6.4 keV line emission observed with 
XMM-Newton (T12). 

3. The spectral analysis of the Arches cluster region, including 
the extended emission around it, reveals a ~1.7 keV 
temperature plasma, which is consistent with the optically 
thin thermal plasma measured by T12 in the cluster’s core 
and associated with thermalization of massive star winds 
and several colliding stellar wind binaries within the cluster 
(T12). 

4. The power-law slope T = 1.62 ± 0.31 and EW 6 . 4 k e v = 
1.1 +q 7 5 keV are in agreement with that measured by T12 
around the Arches cluster and expected for the EECR ion 
model. 

5. The X-ray emission of the Arches cluster is well fitted with 
a model composed of an optically thin thermal plasma and 
a non-thermal component produced by EECR ions (T12). 
The best-fit CR spectral index is s = 1.65+q 5 5 9 . The LECR 
model normalization gives the corresponding power of the 
injected LECR protons in the cloud region at the level of 
(5-8) x 10 38 erg s -1 (with D = 8 kpc and CR minimum 
energy E m { n = 10 MeV nucleon -1 ). 

6. The luminosity of the primary source required by the X-ray 
reflection model is 1-2 orders of magnitude higher than that 


observed in the Arches cluster, which is in agreement with 
similar constraints drawn from the 6.4 keV line flux. 

The X-ray photoionization and CR-induced emission models 
can reproduce the data equally well. Additional NuSTAR ob- 
servations with the Arches cluster on-axis would be useful to 
detect the difference between the emission models above 20 keV. 
At this point, we can conclude that explaining the non-thermal 
emission around the Arches cluster with an X-ray reflection sce- 
nario in which the illuminating source is located in or close to 
the cluster is unlikely. However, irradiation by a long outburst of 
Sgr A* can still be considered, and would imply long-term vari- 
ability of the 6.4 keV line flux should be detected. The excitation 
of the neutral material around the cluster by CR interactions also 
looks plausible, though a large flux of LECR protons is needed 
to produce the observed fluorescent emission. 
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